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Cherenkov Hght from cosmic -ray muons is a significant source of background for the Imaging Atmospheric 
Cherenkov Technique. However, muon events are also valuable as a diagnostic tool because they produce 
distinctive ring images, and the expected amount of Cherenkov light per unit arclength can be accurately 
calculated. We report on a comparison of real and simulated muon events in VERITAS Telescope 1, using this 
comparison to validate the detector model and to determine the light collection efficiency of the telescope. 



1. Introduction 

Determining the absolute energy calibration of a Cherenkov telescope requires knowledge of the signal size 
produced by a single photoelectron. Muons produce sharply defined ring images in the focal plane and provide 
the best "test beam" available for evaluating the end-to-end detector performance! 1 1. Here, we show that 
simulations of muon showers and VERITAS Telescope 1 reproduce the muon energy spectrum, maximum 
Cherenkov angle, and the hght collection efficiency of the telescope. 

VERITAS Telescope 1 f2l is sited at the Whipple Observatory basecamp on Mt. Hopkins, at an altitude of 
1275 m. The telescope has a 12-m diameter reflector and a 499-pixel photomultiplier tube (PMT) camera read 
out by 500-MHz FADCs. The PMTs are Photonis XP2970/02. The pixel spacing is 0.15°. Telescope 1 has 
been taking data since January, 2005|3 1; the 5.6 hours of data used in this analysis were taken during May 9-13, 
2005 at an average zenith angle of ~ 40°. Light cones had not yet been installed when this data was taken. 



2. Simulation and Analysis of Muon Images 

Cherenkov light is produced by charged particles when their speed exceeds the speed of light in the ambient 
medium, f3 > l/n(A). At the altitude of Telescope 1, the threshold energy for a muon to produce Cherenkov 
light is ~ 4.6 GeV, and the Cherenkov angle of a 20-GeV muon is within 0.01° of its maximum value of 1.32°. 

Simulated muons are generated using Corsika 6.03 lj4| for the air showers and GrlSUlSj for the detector 
simulation. Muons are scattered with impact parameters within a radius of 20 m and at angles up to 1° off 
the telescope axis. They are injected at a distance of 2 km (1.5 km vertical distance) from the telescope. 
The atmosphere is modelled using a Corsika external atmosphere appropriate for midlatitude winter (discussed 
below) and atmospheric extinction is modelled using MODTRAN 4|6 1. Measurements of the mirror reflectivity 
are used over the range 260-700 nm, and are extrapolated down to 200 nm. The PMT quantum efficiency (QE) 
is taken from manufacturer specifications covering the full range 200-700 nm. FADC traces which are not 
part of an image are extracted from the data runs and used to provide pedestal traces for the simulations to 
accurately reproduce fluctuations from night-sky background and electronics noise. Individual pixel gains 
and timing offsets are also applied to the simulations. A broken power law is used for the energy spectrum, 
with the spectral indices inferred from Q: indices of -1.93 and -2.34 over the energy ranges 4.5-10 GeV and 
10-20 GeV, respectively. Equal numbers of /i+ and fi^ are simulated. 
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Images are characterized according to a circle fit, minimizing the function 

x^= -(.Vi^^^^ofT^y^^vo? - R)\ (1) 

where Wi is the signal amplitude in pixel i, xq, j/q, and R are the fit parameters of the circle, and Xi and 
are the coordinates of the i*'* pixel in the image. A two-pass fitting procedure is used: after the first pass, 
pixels which are more than 0.3° away from the circle are thrown out, and the image is fit again. This procedure 
improves the fit in cases where a muon is accompanied by light from an associated hadronic shower Figure[2 
shows a muon image recorded by the VERITAS Telescope 1 camera. The black circle is a fit to the image. 

Muon images are selected from the data files on the 
basis of a number of cuts: 

• Hillas width|8| of the image (length of the 
semi-minor axis) > 0.5°, 

• image size (sum of FADC signals in pixels 
included in the 2nd-pass fit) > 3000 digital 
counts (dc), 

• RMS timespread in the arrival of light in the 
pixels < 3.0 ns, 

• RMS width of the Cherenkov ring < 0.12°, 

• completeness of the ring (> 75 % of azimuthal 
bins with a size above a 50-dc threshold), 

• ring center within 0.4° of the center of the field 
of view, and 

• a negligible change in ring parameters be- 
tween the first- and second-pass fits. 

These cuts select complete and nearly complete rings 
with very little contamination from hadronic show- 
ers. Except where otherwise noted, all results below 
utilize all of the above cuts. 

Following the application of cuts, the simulated and data image size distributions are compared, and a scale 
factor fs — <^'^^>"'" is determined. If necessary, a second set of simulated muons is then generated with the 

^ S12^C ^ data 

PMT QE scaled by fg, independent of wavelength. This process is repeated until /g converges to unity. The 
net required change in QE is then a measure of our understanding of the telescope's light collection efficiency. 

3. Results and Discussion 

The results shown below required no iteration. Figure |2] shows that the simulated muons reproduce basic 
characteristics of the data reasonably well. In each case, the simulated curves have been rescaled so that the 
area of the simulated histogram (dotted line) equals that of the data histogram (solid line). In addition, the 



Run: 1252 Event: 13979 Type: 1 GPS: 2005 132 : 9 : 10 : 8.75621 
Max channel 500 



Num Samples 24 
Num Trigger 40 
Num Tubes 65 
Num Dead 28 
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Figure 1. Typical muon image. 
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simulated image sizes (Figure|2^) have been rescaled by fg = 1.0034 so that their mean equals that of the data. 
The distribution of rescaled simulation image sizes agrees well with the data. The distributions of muon angles 
with respect to the telescope optic axis (Figure|2j), with cut at 0.4° removed) and the azimuthal distributions of 
light (Figure|2j;) also agree well, showing that the simulations sample the same distribution of off-axis angles 
and impact parameters as the data. 

The shape of the Cherenkov angle distribution (Figure depends on the muon energy spectrum, and it is 
likely that the muon spectrum is softer than assumed at a zenith angle of 40°. The maximum Cherenkov 
angle depends on the index of refraction near the observation level; the good agreement (within ^ 0.01°) 
between simulations and data on the maximum Cherenkov angle (^ 1 . 32°) is achieved by using an atmospheric 
model appropriate for midlatitude winter A more appropriate model, the US 1976 Standard Atmosphere, 
underestimated the maximum Cherenkov angle by ^ 0.04°. Several possible explanations for this discrepancy 
are under investigation. 




dataSize 


Entries 


1110 


Mean 


3580 


RMS 


360.7 


simJ 


Size 


Entries 


606 


Mean 


3583 


RMS 


345.3 



GOOO 6500 7000 

Image Size (dc) 




Off-Axis Angie (deg) 




Number of Plii Bins Contributing 



(d) 



50 — 
40 |- 
30 — 
20 — 




dataCherAngie 



Entries 1110 
Mean 1.213 
RMS 0.08426 



simCher Angle 



Entries 606 
Mean 1.214 
RMS 0.0799 



1.5 



Clierenkov Angle 



Figure 2. Properties of muon rings. Solid line: data. Dotted line: simulation. Simulated distributions have been scaled to 
have the same area as the data distributions, (a) Size spectrum, (b) Angle between muon direction and telescope axis, (c) 
Azimuthal completeness of rings, (d) Cherenkov angle. 

The dc/pe ratio is calculated from simulations as in Figure|3]by selecting simulated events that pass the above 
cuts and fitting a line to a plot of the size in digital counts (dc) of the ring versus the number of photoelectrons 
generated by the simulation. The dc/pe ratio is the slope of the line, 5.52 ± 0.09 dc/pe, which is in good 
agreement with the value found via special laser runs, 5.3 ± 10 %syst dc/pe. It also agrees with the expected 
value of 5.57 based on the design parameters of the electronics chain, which verifies that the simulation is 
working as expected. 

If we take the dc/pe ratio as well measured, then we can use the scale factor fs to verify our understanding of 
the telescope's light collection efficiency. The telescope's light collection efficiency is determined by several 
factors: atmospheric extinction and scattering, mirror reflectivity, PMT QE, and the packing fraction of the 
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camera. For this work, it is convenient to scale the PMT QE by fs- By iteratively generating simulated muons, 
each time scaling the PMT QE, we can determine the adjustment to the collection efficiency required to achieve 
/s = 1- We find that the base simulation, modified to include the extrapolation of mirror reflectivity over 
200-260 nm, reproduces the telescope's collection efficiency to within ~ 0.5% on the first iteration. An earlier 
series of iterations which neglected this extrapolation 
required an efficiency adjustment of ^ 13 %, consis- 
tent with the indication from simulations that for lo- 
cal muons approximately 10 % of the detected light 
is below 260 nm. Note that for gamma rays the con- 
tribution below 260 nm to the detected light is negli- 
gible. 
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4. Conclusion 

Comparison of muon rings in VERITAS Telescope 
1 data with simulations demonstrate that the simula- 
tion of the VERITAS Telescopes reproduces well the 
properties of the images. Combined with measure- 
ments of the single-PE amplitude via dedicated laser 
runs, this study shows that the base simulation rea- 
sonably reproduces the light-collection efficiency of 
the telescope. 

A dedicated high-multiplicity trigger is planned 
which will allow the individual telescopes of the 
VERITAS array to acquire muon images concurrent 
with regular data taking for continuous monitoring 
of the telescopes' absolute gain and light collection 
efficiency. 
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Figure 3. Image size versus number of photoelectrons for 
simulated muon events. The slope of the line is the dc/pe 
ratio, a measure of the absolute gain of the telescope. The fit 
parameters define a line y = pQ + pi * x. Error bars have 
been suppressed for clarity. 
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